Role of in-plane dissipation in dynamics of Josephson lattice in high-temperature 

superconductors 
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We calculate the flux-flow resistivity of the Josephson vortex lattice in a layered superconductor 
taking into account both the inter-plane and in-plane dissipation channels. We consider the limiting 
cases of small fields (isolated vortices) and high fields (overlapping vortices). In the case of the 
dominating in-plane dissipation, typical for high-temperature superconductors, the field dependence 
of fiux-flow resistivity is characterized by three distinct regions. As usual, at low fields the flux-flow 
resistivity grows linearly with field. When the Josephson vortices start to overlap the fiux-flow 
resistivity crosses over to the regime of quadratic field dependence. Finally, at very high fields 
the fiux-fiow resistivity saturates at the c-axis quasiparticle resistivity. The intermediate quadratic 
regime indicates dominant role of the in-plane dissipation mechanism. Shape of the field dependence 
of the fiux-fiow resistivity can be used to extract both components of the quasiparticle conductivity. 



Stack of low dissipative Josephson junctions formed bsz. 
the atomic layers of high-temperature superconductorsEJ 
represents nonlinear system with unique dynamic prop- 
erties. Magnetic field applied alona-the layers creates the 
lattice of Josephson vortices (JVs).B Transport properties 
of this state are determined by dynamics of the Josephson 
lattice. Two distinct regimes exist depending on strength 
of applied magnetic field. At low fields the Josephson vor- 
tices are isolated and form a triangular lattice, strongly 
stretched along the direction of layers. An isolated JV 
is characterized by the nonlinear core, the region within 
which the phase difference between the two central lay- 
ers sweeps from to 27r. The core size is given by the 
Josephson length, 7s, where 7 is the anisotropy of the 
London penetration depth and s is the interlayer spac- 
ing. This regime of a dilute lattice is characterized by the 
linear field dependence of flux- flow resistivity pff oc B. 
The linear flux-flow branch in BizSrsCaCuaOg (^fiJCO) 
at small fields has been observed experimentallyflo 

When magnetic field exceeds the crossover field, Bcr — 
4>o/7r7s^, the cores of JVs start to overlap and a dense 
Josephson lattice is formed, in which JVs fill all layers.El 
In this regime the linear field dependence of the flux-flow 
resistivity breaks down. Further field behavior depends 
on the mechanism of dissipation. 

Moving Josephson vortices generate both in-plane and 
inter-plane electric fields, which induce dissipative quasi- 
particle currents. Usually only dissipation due to the tun- 
neling of quasiparticles between the layers is rtakcn into 
account in calculation of the viscosity of JVs.l3 However 
in the high-Tc superconductors the in-plane quasiparticle 
conductivity aab is strongly enhanced in superconducting 
state as compared to the normal con|d|Uctivity due to re- 
duction of phase space for scattering3l3, while the c-axis 
component ctc rapidljz" decreases with temperature in su- 
perconducting state. Below the transition temperature 
the anisotropy of dissipation aab / o'c becomes larger than 
the superconducting anisotropy 7^. This leads to domi- 



nating role of the in-plane dissipation in dynamics of the 
Josephson lattice. 

In this Letter we calculate the field dependence of the 
flux-flow resistivity pff{B) taking into account both the 
in-plane and inter-plane dissipation channels. We sep- 
arately consider the regimes of small and high fields. 
The flux flow resistivity at high fields, taking into ac- 
count the in-plane dissipation, has been studied before 
in Ref. 
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For the case of purely c-axis dissipation linear 
growth of the flux-flow resistivity saturates at the c-axis 
quasiparticle resistivity pc when magnetic field exceeds 
the crossover field B^r- Dominating in-plane dissipation 
leads to qualitative change of the shape of pff{B). In 
this case pff also increases linearly at small fields. The 
slope of this dependence is mainly determined by aab and 
at S w Bcr the resistivity p// is still much smaller than 
Pc- At B ^ Bcr the field dependence of p// crosses over 
to even faster, quadratic, dependence. Only at signifi- 
cantly higher field, B w \/'a^dJ{^y^a\^Bcr, Pff reaches 
Pc- Therefore the field dependence of p// can be used to 
extract both components of the quasiparticle conductiv- 
ity. 

Dynamics of the moving Josephson lattice is governed 
by the coupled Siae-Gordon equations for the interlayer 
phase differences. til The equations taking into account in- 
plane dissipation have been derived in Refs. Con- 
sider a layered superconductor in magnetic field applied 
along the layers (jj direction) and carrying transport cur- 
rent along the c-axis (z direction). We express fields and 
currents in terms of the gauge invariant phase difference 
between the layers 6'„ — (j>n+i — i'n — and the in- 

plane superconducting momentum pn = x'Pn ~ ^-^x- 
The local magnetic field Bn between the layers n and 
n + 1 can be expressed as 
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The components of electric field can be approximately 
represented as 
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These expressions are valid assuming fast equilibration 
of the order parameter inside the layers. More general 
situations have been considered in Refs. p^ , p^ p^. The 
components of electric current, and j^, consist of the 
quasiparticle and superconducting contributions 
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where aab and CTc are the components of the quasiparticle 
conductivity, Aab and Ac are the components of the Lon- 
don penetration depth, and jj = c^o/{8tt^sXI) is the 
Josephson current density. Note that the in-plane cur- 
rent is actually concentrated inside the superconducting 
layers and physically meaningful quantities are the two- 
dimensional current densities Jxn m the layers. To be 
precise, the bulk current in Eq. (^) is defined at discrete 
points Zn = ns as jxizn) = Jxn/s. Using above rela- 
tions we rewrite the z and x components of the Maxwell 
equation 
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In the second equation we replaced dB jdz by the discrete 
derivative (i?„ — B„_i)/s. We also neglected the in-plane 
displacement current dDx/dt, because typical frequencies 
involved in Josephson dynamics are much smaller than 
the in-plane plasma frequency c/Xab- Eqs- (0), (||), and 
give closed system which describes dynamics of the 
phases 0„(a;, t), fields i) and momenta p„(x, t). The 

movinglattice generates both in-plane and c-axis electric 
fields (||) leading to dissipation. The rate of energy dis- 
sipation W per unit volume is given by 



V27rc 





(7) 



For the steady state motion with small velocity v the 
phase differences vary in space and time as 9n{x,t) = 
6'ii\x~vt), where 6*1"'' (x) is the static phase distribution, 
and Eq. (M) can be rewritten as 
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where 
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is the linear viscosity coefficient of the lattice per unit 
volume and (...) means averaging with respect to x and 
n. The flux-flow resistivity pj// is connected with r]j by 
relation pjff = B'^/{c^rij). 

Consider the regime of small fields, B ^ B^r = 
i^ol {n^s^). In this regime the JVs are isolated and dis- 
sipation is concentrated in the vicinity of nonlinear vor- 
tex cores. In this case 77,/ is proportional to the field 
•qj = Brjjy/^o, where rjjy is the viscosity coefficient of 
an individual JV per unit length. The viscosity coeffi- 
cient due to the c-jaxis dissipation has been calculated by 
Clem and Coffeyfl Similar problems of viscous friction 
have been studied for an Abrikosov vortex (see, e.g., r£^ 
view |T6| ) and for a Josephson vortex in a single junction.EJ 

In the vicinity of the vortex cores one can neglect 
screening effects and express the phase differences 9n{x) 
and momenta p„(a:) via the in-plane phases 0„(a;), ^i"'' « 
</)„+! - (j)n, Pn'' ~ Vj^c^n (wc are using the gauge divA = 
0) . Numerically accurate phase distribution (/)„ (x) in the 
vicinity of the vortex core was obtained in Ref. nSl 
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with u = X/7S.E3 We can represent now the viscosity 
coefficient rjj^ as 
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Using the phase distribution (^ we compute Cc ~ 9.0 
and Cab ~ 2.4. Finally, we obtain the following re- 
sult for the flux-flow resistivity at small fields pj// = 
<i>oB/(c2,7j,„) 



PJff 
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The case of dominating c-axis dissipatiou. (aab T^o'c) 
has been considered by Coffey and ClemEl. They obtain 
the coefficient 2.8 instead of 4.4 using the approximate 
phase distribution. 

Now we consider the regime of high fields, B > Bcr- fn 
this regime we can obtain a simple analytical result-using 
expansion with respect to the Josephson current .113 For 
the static lattice in the zero order with respect to jj we 
have Bn{x) — B, Pnix) — 0, and = 2ttsBx/^o + 

im. The first iteration with respect to jj = ^^^"^^2 gives 
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Substituting expressions for 6'„ (a;) and p„ (x) into Eq. (g) 
we obtain 



and 
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for B > B^rB At B ^ Bcr this result approximately 
matches with Eq. (px|). Eqs. (fO) and ( ]l^ ) represent the 
main results of this Letter. We see that the shape of 
pjff{B) is determined by the ratio (Jahlil^'^c)- In the 
high-temperature superconductors typically a ah l'^c~l^ 
near the transition temperature. However the ratio 
(^abl<^c rapidly becomes much larger than 7^ with tem- 
perature decrease because of (i) significant enhancement 
of Gab due to-jthe suppression of in-plane |-scattering of 
quasiparticlesQ and (ii) fast decrease of CTc-^ This means 
that dependence pc oc B^ holds in a wide field range 
Bcr < B < Bcr- The field Bcr is almost temperature 
independent and for optimally doped BSCCO (7 « 500) 
Bcr ~ 0.5T. The field Ba- has strong temperature de- 
pendence via (Jab{T) and (Tc{T). Taking typical values 
for T « 20K, = 2 • f O^^^Qhm-cm)"! and aab ~ 
5 • fO^(Ohm-cm)"\ we obtain estimate B„ w 4T. The 
field dependencies for the cases of dominating c-axis dis- 
sipation and dominating in-plane dissipation are sketched 
in Fig. ||. The shapes of the field dependencies are quali- 
tatively different for these two cases. In particular, they 
have opposite curvatures at i? ^ Bcr- Therefore, the 
shape of pjff{B) can be used to extract both compo- 
nents of the quasiparticle conductivity. At present, there 
are no published data for the field dependence of the 
flux-flow resistivity in a wide field range. Dynamics of 
the Josephson lattice rat high fields has been studied by 
G. Hechtfischer et alzJi From the I-U curves presented 
in this paper one can see that the slope dl/dU for the 
first flux-flow branch indeed has a strong field dependence 
with upward curvature in the field range 2 — 3.5 tesla. 



In conclusion, we calculated the flux-flow resistivity 
of the Josephson lattice at small and high fields and 
demonstrated that strong in-plane dissipation qualita- 
tively modifies its field dependence. 
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FIG. 1. Schematic field dependence of the flux-flow resis- 
tivity for the cases of dominating c-axis dissipation channel 
(a-ab/cc ^7^) and dominating in-plane dissipation channel 
{o'ab/o'c 3> 7^). For strong in-plane dissipation the depen- 
dence pff{B) has pronounced upward curvature at B > Bcr 
and approaches pc at the field much larger than the 
crossover field Bcr- 



